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Abstract 

j^ . We present an analysis of the static properties of heavy baryons at next-to-leading 

Tiij- I order in the perurbative expansion of QCD. We obtain analytical next-to-leading order 

three-loop results for the two-point correlators of baryonic currents with one finite mass 
quark field for a variety of quantum numbers of the baryonic currents. We consider 
both the massless limit and the HQET limit of the correlator as special cases of the 
QQ I general finite mass formula and find agreement with previous results. We present closed 

^D ■ form expressions for the moments of the spectral density. We determine the residues of 

physical baryon states using sum rule techniques. 
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1 Introduction 

Baryons form a rich family of particles which have been experimentally studied with high 
accuracy p]. With the advent of new accelerators and detectors many properties of baryons 
containing a heavy quark have been experimentally measured in recent years |lj . A theoretical 
analysis of these experimental data gives a great deal of information about the structure of 
QCD and the numerical values of its parameters. The hypothetical limit iVc — > oo for the 
number N^ of quark and gluon colours in the symmetry group SU{Nc) was especially successful 
for baryons [2]. The analysis of this limit is a very powerful tool for the investigation of general 
properties of gauge interactions. The information about the spectrum of baryons is contained in 
the correlator of two baryonic currents and the spectral density associated with it. Within the 
operator product expansion and to leading order in perturbative QCD the correlator is given 
by a product of Nc fermionic propagators. The diagrams of this topology have been studied in 
some detail [31 HI [5]. They are rather frequently used in phenomenological applications [6l[71[8]. 
In general, these diagrams represent the leading order of the perturbative expansion for the 
relevant correlation functions. In some cases, though (especially for the gluon current correlators 
[9]), they first appear at next-to- leading order. Complete calculations beyond the leading order 
have not been done for many interesting cases. In this paper we fill this gap. 

We report on the results of calculating the as corrections to the correlator of two baryonic 
currents with one finite mass quark and two massless quarks. We present analytical results and 
discuss the magnitude of the as corrections for the physically interesting cases. The high energy 
(massless quarks) and near-threshold (Heavy Quark Effective Theory, HQET) (TUl HH [T21 [13] 
limits are obtained from our general results as special cases. For these we find agreement with 
previous results in the literature. We present analytical results for the moments of the spectral 
density associated with correlators of baryonic currents. 

We briefly discuss the impact of our new results for the baryonic correlators on the phe- 
nomenology of baryons. However, the main aim of this paper is to present the results of the 
perturbative calculations in some detail and to show how they have been arrived at. The new 
technically demanding feature of our calculation is the presence of a finite-mass quark in the 
correlator which is needed for baryons containing a heavy charm or bottom quark. For the 
6-quark the accuracy of the HQET approximation is rather good. The exact result will help in 
controlling the precision of the approximation. For the c-quark, however, the accuracy of the 
near threshold approximation is insufficient for physical applications and the use of the exact 
formulas is unavoidable. 

The massless case has been known since long ago [71 [8] and serves as a test of the massless 
limit of our results. We mention that the mesonic analogue of our baryonic calculation with 
one finite mass quark and one massless antiquark was completed some time ago [14| and has 
subsequently provided a rich source of inspiration for many applications in meson physics. 

Some of the techniques used in this paper have already been usefully employed in the analysis 
of perturbative corrections to sum rules involving pentaquark states p5] . 



2 Generalities 

In this section we present our choice of interpolating currents for baryons containing one heavy 
quark. We also introduce two-point correlation functions as the principal tool in our analysis of 
the static properties of heavy baryons. Finally, we give an outline of the techniques that were 
used in our calculations. 

2.1 Choice of currents and correlators 

A generic lowest dimensional baryonic current has the form 

J = e'''''iulCTd,)T'^a. (1) 

The current ([T]) refers to a baryon with three valence quarks and no gluonic fields and no 
derivative couplings. \E' is a finite mass quark field with the mass parameter m and u and d 
are massless quark fields. C is the charge conjugation matrix, e"''^ is the totally antisymmetric 
tensor and a, b, c are colour indices of the SU{3) colour group. F and F' stand for Dirac matrices 
or strings of Dirac matrices where possible Lorentz indices on F and F' such as in 7^^ or a'^'^ are 
suppressed. In much the same way we have suppressed a possible Lorentz index on ^^ which 
is needed later on in the discussion of the spin 3/2 field. For F = 1, F' = 75 the interpolating 
current has the quantum numbers of a J^ = 1/2+ baryon. Other baryonic currents with 
any given specified quantum numbers are obtained from the current in Eq. ([1]) by using the 
appropriate Dirac matrices or strings of Dirac matrices. We first consider the simplest case 
and take F = 1, F' = 1 corresponding to an interpolating current with quantum numbers l/2~. 
This allows us to explain our techniques and to demonstrate the idiosyncratic features of the 
calculation. Later in the text we will introduce more general interpolating currents and discuss 
calculational differences in comparison to those in the simplest case. 

The correlator of two baryonic currents can be expanded into a basis of invariant functions. 
The form of this expansion depends on the Dirac matrices employed. In the simplest case 
F = 1, F' = 1 there are only two invariant functions Il'^{q^) and n™'(g^) in the expansion which 
are defined through 

tJ{T{j{x)]{0)})e'''dx = mU^{q')+4U''{q'). (2) 

In the following we shall refer to the two contributions on the r.h.s. of ([2]) as the mass and 
the momentum term, respectively. Note that each of the replacements F — * F75 and F' —>■ F'75 
leads to the change W{q'^) -^ —W{q'^). 

2.2 Basic techniques 

The generic correlation function n(g^) has a dispersion representation 

T-r/ 2\ f^ p(s)ds ^ , . 

\l(q)= \ + subtractions (3j 

Jm2 s — q^ 



through its discontinuity p{s) on the physical cut s > m 
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P(') = i^^{^(' + '^)-^('-'^))- (4) 



The discontinuity can also be written as the imaginary part of the correlation function if the 
phases are properly chosen, 

p{s) = -lmU{s + iO). (5) 
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The expression for the discontinuity or spectral density p{s) is simpler than the expression for 
the correlation function n(g^) itself. The knowledge of p{s) suffices for physical purposes and 
allows one to recover the whole function n(g^) through a one- dimensional integral with a simple 
weight function as given in Eq. ([3]). For this reason we concentrate on calculating the spectral 
density p{s). 

The general strategy is rather straightforward. The main part of the calculation is done 
by using a symbolic manipulation program. One reduces all integrals to some basic master 
integrals and then one puts them together again to get the result for a particular correlator 
with any given quantum numbers. This program has been explicitly realized in our evaluation. 

The topology of the NLO diagrams is such that at least one line connecting the initial and 
final points of the diagram is free. If this line is the massive one, the remaining part of the 
diagram consists of massless lines and can be integrated analytically. Adding the massive line 
leads to a one-dimensional integration which can be done analytically. 

If the massive line is part of the radiative corrections, the basic quantity is a NLO two-point 
correlator of the meson type with one heavy and one light quark. The spectral density for this 
NLO correlator is known to be computable. The fact that we only have one finite mass and a 
special topology of diagrams in the baryon sector therefore makes the analytical computation 
feasible. Note that the computation with two different finite masses can still be done but 
requires a numerical calculation while some limiting cases such as the small mass ratio limit 
can still be done analytically. 

In order to explain our main tools, we consider the correlation function of the baryonic 
current, which, up to NLO, can be written as 

U{q^) = J dm2{q - k)U,{k) (6) 

where 11x^2 (^) are one- and two-line correlators. If the one-line correlator is massive, the massless 
two-line correlator can be explicitly integrated and we are left with an integral in D- dimensional 
space-time given by 

(m2 - P)«(-(g - kf 

This integral can be expressed through hypergeometric functions and is therefore completely 
known. If the one-line correlator is massless, the spectral density reads 



^("'/^) / f^2_l,2\a(_(^_l.)2\f3- (V 



p{s)= ds'p,{s')(l-s'/s). (8) 

It is not difficult to obtain p{s) since P2{s) is known from mesonic type calculations. 

All in all the calculation includes no unknown elements in the sense that all necessary 
blocks (prototypes or masters) are known to be calculable analytically. The main problem is 
the reduction of the initial diagrams to prototypes and the assembly of the final results from 
these building blocks. This has been done using the computer. 

In the next section we present the calculation and results for the lowest spin baryons which 
means that the Dirac matrices in the interpolating currents are just unity or 75. 



3 Lowest spin baryons 

In this section we present the results for the simplest choice of the Dirac structure of interpolat- 
ing baryonic currents. We take F = 1 and T' = 75 which corresponds to an interpolating current 
with quantum numbers J^ = 1/2+. We can in fact omit 75 in the process of the calculation 
because the effect of 75 can later on be easily accounted for by a simple multiplication in Dirac 
space. In this case the basic baryonic current has the form 



J 



^abci 



W^C4)^a. (9) 



For this scalar case with F = F' = 1, the results for the invariant functions n''"(g^) and n'^(g^) 
in Eq. ([2]) have already been presented in Refs. [HI [17]. The invariant function n"(g^) with 
a G {q, m} can be represented compactly via the dispersion relation 



p°'{s)ds 
s — q^ 



W{q^) = / ^ ^-^ (10) 



where p"(s) is the spectral density. All quantities are understood to be appropriately regular- 
ized. Since the spectral density is the main object of interest for phenomenological applications, 
we limit our subsequent discussion to the spectral density 

where /x is the renormalization scale parameter, m is the pole mass of the heavy quark (see e.g. 
Ref. [IB]), and a^ = as{ij). 

3.1 LO analytical results 

The leading order two-loop diagram is shown in Fig. [I](a). Note that this topology coincides 
with what is referred to as sunrise-type diagrams for which a general evaluation method (with 
arbitrary masses) has been developed in Refs. [H |19]. Sunrise-type diagrams can be calcu- 
lated by a variety of methods. In this paper we apply the configuration space technique in a 
straightforward manner. The result reads 

p^[s) = l + Qz-Qz^ -z^ + 62(1 + z)\nz (12) 

pI(^s) = --2z + 2z^ -jz^ -?>z^\az (13) 

with z = m? / s. 

3.2 NLO analytical results 

The contributing three-loop diagrams are shown in Figs. [T](bll) to (c21). They have been 
evaluated using the advanced algebraic methods for multi-loop calculations along the lines 



Figure 1: Two-loop (al) and three-loop (bll-c21) topologies with one external momentum. 
Heavy lines represent the heavy quark and light lines massless quarks. 



decribed in Refs. jl]. The result can be obtained analytically. In the MS-subtraction scheme 
one has [Ml [IT] 

m/ N 665 665 9 o /58 , , 9 58 o\ , , 

pf (s) = 9 + --Z - —-z^ - 9^=^ - — + 42^ - 42^2 - —z"" ln(l - z) 
9 9 V 9 9 / 



+ 2 + 



l^^ _ ^^2 _ 58^3\ 1^ ^ ^ 8 ^^ ^ g^ _ g^2 _ ^3^ ^ i2(z) + - ln(l - z) In z 

+z {2A + 36^ + ^z^^ {hi2{z) - C(2) + ^ In^ z^ + 24^(1 + z) (hh{z) - C(3) - h^i^^z) In z) , 

(14) 



pI{s) 



71 565 7 



625 



109 



48 36 



z z H 2; z — 



1 /49 116 116 



36 



48 



49 



36 V36 9 



z H z z ln(l - z 



9 



36 



/I 17 9 113 o 49 .\, 2 / 3 4\ / . ^ 1 . 

+ f - - y^ - ll;z' + —z"^ - —z^\ In^ + - (1 - 8^ + 8;z^ - /) f Li2(;z) + - ln(l - z)\nz 

Af (54 + 8z - z^) (Li2(^) - C(2) + \ In^ z) - 12^2 (^Lialz) - C(3) - ^Li^lz) ln(z)) (15) 



3 V y ^ -V / .V / 2 

where 2; = rrP' j s and Li„(2;) are polylogarithms 



Li„(2;) = Y, 



00 ^fc 



fc=i 



A;*^ 



00 -| 

Lin(l) = E ^ 
fc=i '^ 



CH 



(16) 



C(n) is Riemann's zeta function. Note that in the physical region we have z < 1. Therefore, no 
analytic continuation is required. Note the repeated appearance of particular combinations of 
polylogarithms and logarithms in the form 



Li3(^)-C(3)--Li2(^)lnz 



1 r^ rl7' 

- ^{2U,{z') + Hl-z')\nz') 

6 Jl z 



1 r^ d 7' 

U^{z) - C{2) + -\n' z = J^—{\nz'-Hl-z')) 



and 



Li2(2;) + -ln(l-2;)lnz 



2 Jo \l-z' z' 



(17) 



appear as a consequence of the integration of the two-line spectral functions. 



4 General baryon case 

The baryonic current defined in Eq. ([1]) has the most general structure concerning possible 
choices for the Dirac matrices F and T'. However, not all of this structure has to be kept in 
order to calculate the different invariant functions of the correlator. In four-dimensional space- 
time where the initial (bare) current is defined, Fierz rearrangement can always be used to take 
the heavy spinor out of the trace in the correlator. When taken out of the trace together with 
the heavy spinor, the matrix T' can be considered to be an overall factor that has no effect on 
the calculation. Indeed, 

J = e«''^(Mrcrde)r'^, = r'jr (is) 

where jr is the current for F' = 1, 

jr = e<^^^(M^CF4)^a. (19) 

The general case is recovered as 

n(g2) = z|(T{j(x)j(0)})e*'^^dx = F'nr(g')f' (20) 

with f = 7°F"'^7° and 

nr(g^) =tJ{Tjr{x)MO))e''''dx. (21) 

Because the general result can easily be recovered, we can limit ourselves to the case F' = 1. 
The calculation of the correlator has to be done with its full dependence on the matrix F. We 
obtain the general expression for the correlator in the form 

Ur{q') = J:Aiq')tT,iT,q') (22) 



1=1 



where 



tri(F,g2) 


= Tr(F^f^)m/g2 


tr2(F,g2) 


= Tr(F7J7")m 


tr3(r,g^) 


= TriV^mi/q' 


tr4(F, g2) 


= Tr(F^f7„)7" 


tr5(F,g2) 


= Tr(F7„f^)7" 


tr6(F,g2) 


= Tr(F7J7°)^. 



(23) 

The trace in fl23|) is to be taken only with respect to the 7-string in the round brackets. We 
have again omitted possible Lorentz indices on F such as appear in the next example. 

As an example, let us exhibit the structure of the expressions for the interesting and im- 
portant case of the "vector" current 

f = e'''''iulCYd,)^a. (24) 

In the above notation this means that F' = 1. The expansion of the correlator reads 

m-^iq^) = % f{Tj^'{x)fiO))e"'''dx. (25) 



The correlator can be expanded along a set of ten covariants. The expansion reads 

+i {A\q^y + Alq^g^^ + AUq^Y + AU^^^ + Alq^^^^) / q' ■ (26) 

All invariant amplitudes A" have been calculated. The above expansion of the vector correlator 
is only one of several possible expansions. By making use of the symmetry properties of the 
result we shall later on use a different expansion in terms of nine covariants whose mixing 
property is simpler. 

4.1 LO analytical results 

Differing from our earlier calculation [161 ^\i where the mass and momentum parts were cal- 
culated separately, we have now learned to calculate the mass and momentum parts in one go. 
The reason is that no explicit traces have to be taken when extracting the mass or momentum 
part. Instead, the traces are kept to the very end. It is then not difficult to interprete one part 
of the expression as the mass part and the other part as the momentum part depending on the 
occurence or absence of an explicit factor of m. The result for the leading order diagram (al) 
in D = 4 — 25 space-time dimensions reads 

Pai{z) = EPaiWtr.(r,.). (27) 

Here ^(1, 1) = G/e and G = T{1 + e)T{l - ef/T{2 - 2e). Note that we use hatted spectral 
functions whenever we present them as a function oi z = m? / s. The corrections are 

pI,{z) = D|p^(l,£-2;^) + (l-2;)2pv(l,£;^)-2(^^^^-2;)/)v^(l,£-l;z)|, 

PaiW = -|py(l,£-2;2:) + (l-2;)2p^(l,£;2;) + 2(^^^^ + 2;)py(l,£-l;2;)|, 

PaiW = ^^^Yv{l,e-?>-z) + (^^^^ + z^{l-zfpv{l,e-z) + 

/6-D \ , , /2-D 2Dz ^\ , A 

PaiW = --\pv{\e-?,]z)-{l-zfpv{l,e]z) + 

+ (1 + 3z)/)y(l, e-2-z)-{l-z){l + ?,z)pv{l, e-l-,z)V 

pIM = -Upvil,e-3;z)-{l-z)'pvil,e;z) + 

+ (1 + 3z)pv{l, e-2;z)-{l-z){l + 3z)pv{'^,e -l;z)], 



pIi{z) = -]^pv{l.e-?,-z) + {l-zf{l + z)pv{l,e-z) + 

flD - 10 \ , 

+ ( ^_^ + 32:j py (1, e-2-z) + 

7D-10 3L'-4 o\ , A 

-p3^ + 2-^^^^ + 3^^j/)v(l,5-l;z)|. (28) 

For later reference we also need starred elements p]^i{z), p^ai{z), p^ai{z), p^i{z), p^ai{z), and 
p^ai{z)- These can be obtained from the corresponding unstarred elements by replacing e in the 
argument of the spectral functions pv by 2e. The basic spectral functions pv are given by 

pv{n^,n,-z) = ,rnin/9 T /V - ^)''/'~"^"V-i(x - ;.)^/2-"^c?a: 

r(ni)r(n2)r(l + D/2 -ni- n2) Jz 

(29) 
where T{z) is Euler's gamma function. Note that the singularity given by the factor G(l, 1) in 
Eq. 0271) cancels against the singularity of the second gamma functions in the denominator if 
the appropriate arguments n2 = s — n (or n2 = 2e — n in case of the starred elements) with 
n > are used. In the limit D = 4 we find 

Palis) = ^,T. P:\imy s)tT,{T,s) (30) 



where 



Pliiz) = l + lz-3z'' + z^-^z^ + 2z\nz 

p^l^Jz) = -H z z^ ^H z"^ + (z + -z'^]lnz 

HaiK J g-Tg 2 224V2' 



pIUz) - 


1 

10 


1 ^ , ^2 -,3 I 1 ^4 Is 

2'^' ' ^2' 10^ 




pIUz) - 


1 

40 


1 1 2 1 3 1 4 Is 

a' g' ^2' S' ^QO' 


1 2, 

-z Inz 


pIUz) - 


1 
40 


1 1 2 1 3 1 4 1 5 
4' 6' ^2' S' +60^ 


1 2, 

-z \yiz 


C{z) - 


1 
40 


1 1 2 1 3 1 4 Is 

z z + z z + z 
4 6 ^2 8 ^60 


l.M„. 



(31) 

Note that p'^{z) = Pai{z) = Pai{z). Eqs. ( l30l) and ( 1311) give the full answer for the leading 
order contribution to the baryonic correlators for all possible configurations of Dirac gamma 
matrices. In this sense this completes the leading order calculation of correlators for baryons 
with any quantum numbers as long as there are no derivative couplings in the interpolating 
currents. 

4.2 NLO contributions 

The NLO contributions result from four different diagrams. In the calculation of these diagrams 
we have used different techniques depending on their topologies and on the location of the 
massive line (cf. Fig. [1]). 
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• The NLO light contributions result from two diagrams, the self energy correction of one of 
the massless lines (b21) and the diagram with gluon exchange between two massless lines 
(ell) which we term "light fish". The technique of calculating these diagrams consists in 
first analytically calculating the massless part and then adding the massive fermion line. 
Advantage is taken of the fact that massless two-loop diagrams are explicitly calculable. 
In this case one cannot avail of a convolution of spectral functions. 

• The massive line self energy contribution (bll) was calculated again by using the explicit 
evaluation of the massless part and the dispersion relation for the mass operator of the 
massive quark at leading order. Here we define the heavy quark mass through its pole 
mass which is most convenient for the calculation using cuts. 

• The most demanding semi-massive fish diagram (c21) was calculated by making full use 
of the decomposition into prototypes and the convolution of spectral functions. For this 
diagram the use of symbolic manipulation programs is indispensable as the number of 
terms in intermediate expressions involving different structures are quite large. 

These three main parts constitute the whole calculation. 

4.3 Light contributions (b21 and ell) 

The NLO light contributions result from two diagrams, the self energy correction of one of 
the massless lines (b21) and the diagram with gluon exchange between the two massless lines 
(ell), called "light fish". More details on the calculation of these two contributions are found 
in Appendix A. It turns out that the dominant singular part is proportional to the leading 
order contribution. The result for the spectral density reads 

piight(s) = — — — \['^ + y) ^"^^^^ [~ I I 

where 

(^2 3D/2-4 

Pai(^) = 2{L)^D {Paii^Vs)tT,{r,s)+plUmVs)tT,ir,s)) (33) 

and where Pai{z) are the starred elements introduced after Eq. fl28l) . The spectral functions 
pZ*iz) = -py(l,2e-l;z), 
pIUz) = -^{(l + ^)py(l,2£-l;;.)+py(l,2£-2;z)} (34) 

are known from the scalar calculation, and 



50 Bi / 1 1 17\ Cb 



cf / 1 1 1\ /I 19 






^fi + ^Ui 

6 \e 2j 



(35) 
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where Cp = {N^ - l)/2N^ = 4/3, Cb = {N^ + 1)/2N^ = 2/3 for A^^ = 3 colours. The occurence 
of Cp in the denominator of the second parts in Eq. (155]1 results from the fact that Cp is 
factored out in Eq. (15^ . The parameter cr is defined by 

7af7" = erf, cr = sr{D - 2rr) (36) 

and can be expressed by the signature sr = ±1 of the matrix according to whether one has an 
even or odd number of 7-matrices (including 75) and the number rp of Dirac matrices other 
than 75. The fact that the most singular term is proportional to the leading order contribution 
allows one to extract a common renormalization factor. The result reads 



18(47r 
c 
Air 



f^ (C'^(^)#(mV.) + C^^'^giWmVs) + C-(°)#)(mV.)) [ (37) 



where 



C"^(o) = 2BoG, C"'^^) = 2 I Bo In ( ^ j +Bo + Bi\g, C"'(°) = 2B[G. 
The definition of the functions g'^{z) can be found in Appendix A. 

4.4 The massive line self energy contribution (bll) 

We divide the leading order contribution in Eq. (1271) into a mass part 



(3^ 



{D-2)G{l,l)mn-2,m,. 



pZis) = %;^)%yZTf^"-^/^-(-V^), where 
PTiim'/s) = J:pUmys)tniT,s), (39) 



i=l 



and a momentum part 

nUs) - (-P-2)G(l,l)Arj ^_2^ 2/ X 1 

Pai[s) - ig(47r)i5(/) _ 1)2 ^ Pai[m/s), wucre 

pliimVs) = J:pU^ys)tniT,s). (40) 

i=3 

The result for the self energy correction of the massive line is obtained by a convolution of these 
leading order contributions. One obtains pbii(s) = Pbii(s) + plnis) where 

P'''^'> - 16(4vr)^(D - 1)2 A [~^(r^^^^'^''>di)P'^'^'/''''^''' 

P'''^'' - 16(47r)^(D - 1)2 A 1^ x{l-x)^^^'^'''dx)^-'^'^''''^'' ^^^' 
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and where again z = m? / s. We have introduced the functions 

'D + 2 D-2 



Pa{z) = (l + z)py(l, l;z), 

[pa+b{z) = Pa{z) + pb{z)y and 



Pb[Z) 



z) pv{l,l;z), 



LJz) 



Pb{z')dz' 



(42) 



(43) 



/o {l-z'y 
The details of the calculation can be found in Appendix B. 

4.5 The semi-massive fish (c21) 

The two diagrams obtained by connecting the massive line with one of the massless lines 
(see Fig. [T](c21))are called semi-massive fish diagrams. Summing up the results for these two 
diagrams we obtain 



Pc2l[S) 



g2,D/2-3 



16(47r)3^/2(D-2)(D-l) 



dsis^ {s - Si)py{1,1;si/s) X 



X 



4{D-2){D-2 



Dsi 



(pUzi) tri(r, Si) + p'^izi) tr;(r, Si)) + 



+4(D - 2) (^1 - -jj {pm{z,) tr2(r, Si) + p'^{z,) tr'2(r, s,)) + 
+ (D-2 + 2{D-2f-^ + {D + 2fJ\ (p,(zi) tr3(r, si) + p'^{zi) ti'.iT, s,)) + 
+ j^ (p,(zi)- (tr4(r, si) + tr5(r, si)) + p'^izi) tr^F, Si)) + 
Pg(^i)2 (^^4(r, si) + tr5(r, si)) + Pg{zi) tr'r^{T, Si)j + 



+ 1 



Sl 

s 



- 1- 



+ 1- 



s 

Si 



1 + 



Si 



Pq{zi) tTe{T, si) + pJzi) ti'f^iT, si)) + 



+4P _ 1)^ (^p-(^^)l (tr4(r, si) + tr5(r, si)) + p'^'iz,) tr^(r, si; 



Pm(^l) 



52, 



where zi = m? / si and 

p^{zi) = 2{D-l)B,, 
p,{zi) = 2{D-2){{D-2 + Dzi)B4-2{D-l)ziB6)-D{DB5 + {D-2)Br) 
p;'(zi) = 2{D-2){{l-z,)B, + {D-l){Bs + z,B,)) + D{B, + {D-2)Br), 
p'^izi) = DB, + {D-2)Bj, p'^'iz,) = -B, 
with 

B, = -pv{l, 0, 0, 1, 1; z,) + 2(1 - z^)pv{l, 0, 1, 1, 1; Zi) + 



(44) 



(45) 
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+(1 + z^)pv{l, 1, 1, 0, 1; z^) - pv(l, 1, 1, 1, 0; z{) + (1 - z^fpv{l, 1, 1, 1, 1; ^i), 
B2 = -pv{l, 0, 0, 1, 1; z^) + py(l, 0, 1, 0, 1; z^) - py(l, 0, 1, 1, 0; z^) + 

-py(l, 1, 1, -1, 1; z^) + pvil, 1, 1, 0, 0; z^) - (1 - z^)pv{l, 1, 1, 1, 0; 2i), 
54 = -2/3y(l, 0, 0, 1, 1; zi) + 2(1 - z^)pv{l, 0, 1, 1, 1; zi) + 

+2(1 + z{)pv{l, 1, 1, 0, 1; 2:1) - 2py(l, 1, 1, 1, 0; Zi) + (1 - Zifpv{l, 1, 1, 1, 1; zi), 
B, = 2py(l, -1, 0, 1, 1; ^1) + 2py(l, 0, -1, 1, 1; z^) + 

-2pv(l, 0, 0, 1, 0; zi) + 2(1 - Zi)py(l, 0, 0, 1, 1) - 2(1 - z,)pv{l, 0, 1, 0, 1; z,) + 

+2(1 - ^i)py(l, 0, 1, 1, 0; zi) - 2^ipv(l, 1, 1, -1, 1; z^) + 2(1 + zi)pv(l, 1, 1, 0, 0; zi) + 

-py(l, 1,1,1, -1; z,) + (1 - zi) V(l, 1, 1, 1, 0; zi), 
Se = 2pv(l, 1, 1, 0, 1; z,) - py(l, 1, 1, 1, 0; z,), 
Bj = -2pv(l, 0, 0, 1, 1) + 2pv(l, 0, 1, 0, 1; z,) - 2/3^(1, 0, 1, 1, 0; z,) + 

-2pv(l, 1, 1, 0, 0; zi) + py(l, 1, 1, 1, -1; z,) - 2(1 - Zi)pv{l, 1, 1, 1, 0; z^), 
Bs = 2pv{l, 0, 0, 1, 1; zi) - 2z^pv{l, 1, 1, 0, 1; zi) + pv(l, 1, 1, 1, 0; zi). (46) 

The spectral functions pv'(^i,''T'2,^3,^4,'^5; -2^1) are so-called prototypes which are defined in 
Appendix C The additional traces that do not appear in Eq. (123|) are given by 

tr;(r,g2) = ^^{Tv{Ti[a,^J]i)a^^ -a^^i:v{[a,^J]m)}, 
tr^(r,g2) = ^{Tr(r7„[a^.,r]7")a'^''-a'^'^Tr([a^„r]7.r7")}, 
tr'3(r,g2) = ^^{TT{Ti[a,.,T]i)ia^^-a^''Ti{[a,^J]iTi)i], 

tr;(r,g2) = l{Tr(r^a,,f7j7"^'^^-7"Tr(r7.ra^.^)a'^^ 

+a'^^Tr(ra^,^f7,)7° - (T^'^7"Tr(r7„f^a^,)}, 

K(r,g') = ^{7"Tr(r7„a^.r^)a'^^-Tr(r^fa^,7a)7V^'^ 

+^^^^7" Tr(ra^,7,f ^) - a'^^ Tr(r^f 7,a,,,)7"}, 

tr;,(r,g2) = i{Tr(r7a[a,„f]7")^a'^'^-a'^^Tr([a,„r]7.f7")4. (47) 

The traces tr^(r,g^) correspond to traces which are absent in the leading order term. They 
enter the calculation in the course of renormalizing the results. 

4.6 Renormalization 

The baryonic currents need to be renormalized. In general, there will be mixing under renor- 
malization and therefore one has to construct the whole matrix of renormalization constants. 
Within our technique of parameterizing the results for arbitrary F matrices this is a straight- 
forward procedure. The genuine first-order vertex correction for the current j is given by 

3^ = e'^'^u^^Ci^.^^T + r7.7^)4)7^'^r'vl/, (48) 
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which (using 7/^7^ = gf^" — ia^'^) can be written as 

= e'''"'{2D{ulCTdt)T'^, - {u^C{a^,T - ror^,)4)a^"r'^,}. (49) 

In calculating the correlator we therefore expect objects of the form 



AD- 



:{J^J^) 



Tr(r o r 



'2D ^^^'^ ° ^"^^"^ " ^""''"^ °^ ° '"''' + 



+^^^" Tr((or,,,r - Ta^,) o f o) o . 



(50) 



The open circles stand for further Dirac structures in the calculation. For example, if the first 
term on the r.h.s. of ( l50l) is given by Tr(r o F o) o = g^ tra = Tr(F^r^)^, the whole right hand 
side of Eq. ( 150|) reads 



Tr 



Tr 

2D 



o-fMi^T - Tafj.^] 



ia 



flU 



+ ^cr^''TT[ia,.r-Ta,. 



^^tra-^tr;, 



(51) 



where we have used the trace definitions fl23l) and P7|) . The left hand side of Eq. fl50|) represents 
the singular contribution of the diagram. On the other hand, the singular parts of the spectral 
functions of the basic structures tr^ are 2as/37Te times the LO term, whereas the spectral 
functions of the primed structures tr^ are —ag/Qire times the LO result of the corresponding 
basic structure. Note, finally, that we need to consider only the LO e singularities within the 
MS-scheme. Therefore, we need not specify F at this point. If we write the total result for the 
semi-massive fish as 



E[p^tr,+prtr: 



6 



E 



Poo + Poi^ + -(p1o7 + Pli)) tr,, +- (Pl^; + p1^) tr^ 
7f ^ e 'J n ^ e ' 



(52) 



with p\q = |pqq and p*/q = — |poo5 "we can extract the renormalization factors and obtain 



6 r 

E 



1 + ^) (poo + Poi^ + V (/'n - f/'oi) ) tr. + 



vr 



6,-Poo + ^Pr.)tr: 



(53) 



p; 



Poo + Poi^ + — (Pii-^Poi 



(54) 



Up to 0[oL^ we have 

■(Pii- ^, 

If we substitute this expression for Pqq in the coefficient of tr^, the next-to-leading order con- 
tribution can be skipped while the term proportional to e leads to a subtraction of the finite 
term. 



6 

E 

1=1 



1 



2a, 



p\ tri + 



bixe IT 



1 



^{Pn + -QPo.]tr[ 



In total we have 



Pi tr. +pr tr: = (1 + ^) Pi tr, + (-^p^ + p^'j tr^ 



(55) 



(56) 
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where 



^(A + k)- 



pt^-{Pu + M- (57) 

The coefficient of the primed structure has no LO contribution. Nevertheless, the renormahza- 
tion works in the same formal manner as for the basic structure. As these calculations show, 
the renormalization of the mixing operators for the basic and primed structure is accomplished 

by a renormalization matrix, 



Pb 
.Pt. 




(58) 



Q-KS 



One can easily invert this renormalization matrix to compute the renormalized quantities. One 
also needs the bare quanties which are given by 



Pr 



1- 



2a. 
Svre 



Pi 



P'r 



pi; + 



a 



■5 J 



67re 



Pb- 



(59) 



The results given in Appendices A to C are already renormalized ones. In order to obtain the 
total result we have to combine the spectral functions giving 



Po(^) 



Pleading l^j 



P\\^) = Plight W+P: 



' iW+PfishW, 



Pl'(^) 



Pfish(^)- 



(60) 



Using these spectral functions and the traces defined in Eqs. (I23l) and (1471) . the spectral density 
is given by 



P{s) 



5127r4 



E 

1=1 



Po(^Vs) 1 + 



' ^ 3-4rr + r^\ a« a. 

^m + 5 — tli - — 

6 it: DVT 



tr' 



In 



m 



P^ 



+ ^(/)i(mVs)tri+/31'(mVs)tr:,) 



+ 



(61) 



where n^ = 0, 1 depending on whether there is a factor of m in tr^ or not. Explicitly one has 

3-6 - 0. 



n]^ = 1 and n^ 



5 Some properties of the spectral densities at NLO 



In this section we consider two limiting cases of the NLO spectral densities. First we consider 
the large energy or equivalently the mass zero limit. Second we analyze the near threshold 
limit relevant for a comparison with HQET results. Both limits are interesting and physically 
relevant. The limiting cases for the lowest spin baryons have been discussed before in Refs. [121 
[T7] . We therefore concentrate on the case of the "vector" current (1211) in the following. 



5.1 High energy expansion 

In the high energy (or, equivalently, small mass) limit z —>■ the spectral density reads 



P'^is) 



5127r4 



a.« 



M^ + ?ll2+'"lT 



59 , /p^\\l q'^q" of. "s /13 2, //i2 



'"^-^^ + f T + 3'" T '"»' + 



16 



4 I , ds 

5 [ TT 

if as 

5 TT 



343 1, 

m 

180 3 



s 



313 _1 //i^^ 
180 3 ^ \ s , 



4 



a 



TT V 60 V ^ 



82 



1 



7V-7 1 + — 77-77ln — 



vr 



TT V36 3 \ s 



2m^ 



7^g'' 



^^7*^ 



fY 



45 

vr 



3 

'9 



+ ln 



/^^ 






m^y'^Y 



(62) 



where ?ti = Tny^[n) is the MS mass. When one compares ( l62l) with the results of an afe initio 
(multiphcatively renormahzed) massless calculation one does not obtain full agreement. We 
discuss an alternative route. Instead of expanding p^^'^ along the set of ten covariants given in 
Eq. fl26l) we exploit the symmetry of the problem and expand along an alternative set of only 
nine covariants 



In this case we obtain 

„2 



p^'^{s) 



5127r4 



A. as 59 , 

4 n + — — + In 

I TT I 12 ^ 



4 I as 

5 [ 71 

if, "s 

5 TT 



343 

^180 

'313 

180 



3 \s, 

3 U . 



/^^ 



771- 



qf'q'' 



H-vM'" 



,gV 



^Mg- 



if a» 

5 TT 



313 1, 

m 

180 3 






5 1 TT 1 45 3 Is, 



TT \36 3 Is 



m- 



fYi 



J 



i^ji^y 



a.. 



vr 



9 3 I s 



- m(?^^ + 



g'^7*' + 



/ii/ 



+ 



771'-)^'-)'^ 



(63) 



where the first six contributions which contain both LO and NLO parts are reproduced by the 
massless calculation. It cannot be expected that a multiphcatively renormahzed massless result 
can reproduce the remaining three contributions. In addition, the massless calculation cannot 
reproduce terms such as 2;ln(z) as they appear for instance in the full NLO results in Eqs. ^\M 
and ( TT5I) . These terms can be parametrized with condensates of local operator^. 

5.2 Near- threshold expansion 

As concerns the A-type baryons we have already compared the near-threshold limit with the 
HQET result in Ref. [16]. Here we repeat the exercise for the S-type J^ = 1/2"*" baryon with 
r = •y^, T' = 7^75, i.e. sp = —1, and rr = 1. In the near-threshold limit -E — > (s = {771 + EY) 
Eq. (1^ gives 



Pthrim,E) 



E' 



AOtt'^771 



as /42 4:71^ 

2U + — — + + 2h, 

1 TT I 5 9 \2E 



as fl21 47r2 

1 + — + + 2n, 

TT \ 15 9 \2E 



fJL 
f_p\ 1 



(771 + 



q^'q" 



+ 



+ -In 
3 



771 



P^ 



[m + i)g^'' + 



^For a discussion of this point cf. Ref. [IB] 
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611 OTT 



-S— "(^ 



1 -In 



«,s 



m 



f^' 



m- 



fYi 



+ 



m 



where m is the pole mass. The HQET contribution is obtained by projecting on both sides 
with (1 + t^)/2 = (ttt, + (ji)/2m. Omitting the projector (1 + l>)/2 itself we obtain 



Phqet("^,^) 



E^ 



2QTX^m 



TT I 5 9 



+ 21n(-^ 
\2E 



as 



a^(f(f_ 
Stt s 



1-ln 



2E 



s 

2\ ^ 



+ 



21„|AUilnf!^ 



m 






K^. 






7^7 + 






g^^ + 



qf^Y 



m 



+ 



YY + 



m 



(65) 



We now proceed to compare (p5|) with the corresponding result derived from the HQET cur- 
rent [22] 

hi = iq^Cylqhil5Q, (66) 

where 'jj_ = '^uid^'^ ~ q^q'^ /q^)- One needs to extract the transverse piece from Eq. (165!) using 
the transverse projection g^^ — q^qv/q^- 



mpY,i{rn,E) 



3E^ 



1 + 



tto 



136 

"15" 



2Qn^ 1^ TT 
This can be compared with the HQET result 

3E^ 



pME^fJ') 



2077^ 



1 + — 
TT 



47r2 
9 

116 



2 In 



An' 



2E 



3 



m 



fi^ 



9 



21nf^ 
\2E 



(67) 



(6^ 



where the explicit mass factor m appearing on the left hand side of Eq. (1671) has been absorbed 
in the definition of p-£i{E, jj) in Eq. (l68l) rl After addition of the matching coefficient [22] 

Gsi(m//i, «,) = 1 + ^ f 2 - In l"^] ] (69) 



/^ 



one obtains mp-si{m, E) = Csi{rn/p, as)^psi(yU, E). In this case the matching procedure allows 
one to recover the near-threshold limit of the full correlator starting from the simpler effective 
theory near threshold [21]. Note that the higher order E/m corrections to Eq. (1M|) can be 
easily obtained from the explicit result given in Eq. (IT^ . Indeed, the next-to- leading order 
correction in the low energy threshold expansion reads 



Apj:i{m,E) 



UE^ 



120m2 



as fUM 477 
1 + — Kt:^ + 



120m2 



1 + 



n \ 165 
a, /'512 



n 



51 



9 

47r2 
~9~ 



-^-'^MKi) 



3 

28 



p- 



m + 



In 



m 



51 V /W^ / 
^Note that we have to take into account the contraction with F' ■ • • F' as weU. 



110, / m 

H In 

51 \2E 



I (70) 
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To obtain this result starting from HQET is a more difficult task requiring the analysis of 
contributions from higher dimension operators. 

5.3 Interpolation 

We now discuss some quantitative features of the correction given in Eq. fl6Tl) . Of interest 
is whether the two limiting expressions (the massless and threshold limit) can be used to 
characterize the full function for all energies. To this end we compare components of the 
baryonic spectral density up to NLO. In Fig. [2] we show the ratio p\{s)/ Pq{s) for i = 3. Fig. [2] 
shows that one would obtain a rather good approximation for the full NLO order correction for 
all values of s if one were to use an interpolation between the two limiting cases. 



,25 



if) 

ro o 



^20 
<f) 
a 

10 
5 


-5 



I ' ' ' I ' ' ' I ' 
exact solution 
threshold limit 
massless limit 







2 



4 



12 14 16 
s (in units of m^) 



Figure 2: The ratio p\{s) j p^{s) up to NLO as a function of the squared energy s. 



5.4 Moments 

An instructive set of observables are the negative moments of the spectral density (see e.g. |25j rl. 

'2n 6 



M-n 



s "■p{s)ds 



m 



5127r4(n - 3) fr^ 



5:(M;_,tr,+M:'_,tr' 



-4 ^'■t T^J''-'„_4 ^'■i ; 



Mf) = (n + 1) / z''p''^'\z)dz 
Jo 



(71) 



where p{s) is taken from Eq. (|6T1) and p\rn? / s) and p^'{m^/s) are the coefficients of trj and tr^. 
One has 



Ml 



Mf^ 1 + ^ Ul + 



3-4rr + rr\ ( p^ 



TC 






^We have normalized M^ such that Af^ — 1 for p^{z) ~ I. 



19 



Mil = M;(°)^-^lnf^')+^^nK (72) 



where 



Mm 

n 

Mm 



Gvr V uv' / vr 



12 



(n + 2)2(n + 3)(n + 4)(n + 5)' 

6 
(n + 2)2(n + 3)2(n + 4)(n + 5)' 

12 
(n + 2)(n + 3)(n + 4)(n + 5)(n + 6) ' 



M4(o) = ^ = M^(°) = M6(°) r73l 

(n + 2)(n + 3)2(n + 4)(n + 5)(n + 6) « « ^ ^ 

and (5; = (J^^'') + Sf^ + (5;(^) with 

53(a) ^ 24+ 2N (^ , . + 1) + 7p ''^^'' + ^^ ^ + -rr = (J^^'^) = ^^^'^^ = (5''('^\ 

3^ ^ ^^1^^^ ^ ^"^ 4(n + l)(n + 2)y 3^ " " " ' 

3n3 + 8n2 - 3n - 12 ,,, ,\ 47r2 ,. 

2(n + l)^(n + 2)(n + 3) / 9 

,0. ^ 4 /n^ + 13n^ + hhn^ + TSn^ _ 26n - 82 , , , 

^n = 7T L ^.w ^^w ^^w ^,w ^.^ OA^+1 +7i^ + 

3 \ (n + l)(n + 2)(n + 3)(n + 4)(n + 5) 

n^ + IGn^ + 122n^ + 448n3 + 655^2 + 120n - 300 ,,, A An"^ .,. 



4(n+l)2(n + 2)(n + 3)(n + 4)(n + 5) ^ "■ ') 9 

^4(,) 4/ -2(n3 + 9n2 + 24n+19) ,,, ., , 

3 y (n + 1) (n + 2) (n + 3) (n + 4) 

n^ + 8n5 + 16n^ - 12n3 - 21^^ + 108n + 132 , ,, ^A An"^ .,. 

1! (n -I- 2 I -I = ^^ ' 

4(n + l)2(n + 2)(n + 3)(n + 4) ^^ ^J 9 

n2-2n-2' 



^(n+l)+7i;+ ,, , ,, = C 



3(n + 4) \^--v- ^ ■ '- ■ 4(n + l 

Further ip{z) = r'{z)/r{z) is the digamma function, ip'{z) its first derivative (first polygamma 
function) and 'Je = 0.577. . . is the Euler-Mascheroni constant. In order to combine the three 
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n 


rr = 


(scalar) 


rr = l 


(vector) 


rr = 2 


(tensor) 




Al 


S'n - SLl 


Al 


51 - SLi 


Al 


€ - sLi 





2.66667 


— 


3.33333 


— 


4.00000 


— 


1 


5.59259 


2.92593 


4.92593 


1.59259 


5.14815 


1.14815 


2 


7.39815 


1.80556 


5.98148 


1.05556 


5.95370 


0.80556 


3 


8.70741 


1.30926 


6.75741 


0.77593 


6.55185 


0.59815 


4 


9.73056 


1.02315 


7.36389 


0.60648 


7.01944 


0.46759 


5 


10.5674 


0.83683 


7.85786 


0.49397 


7.39913 


0.37968 


6 


11.2735 


0.70611 


8.27230 


0.41444 


7.71635 


0.31722 


7 


11.8831 


0.60957 


8.62791 


0.35561 


7.98730 


0.27095 


8 


12.4186 


0.53552 


8.93842 


0.31052 


8.22281 


0.23552 


9 


12.8956 


0.47701 


9.21341 


0.27499 


8.43046 


0.20765 



Table 1: Values for the rational part A^ of the first moments 5;^ and their relative difference 



^n ~ ^n-i fo^ ^he scalar, vector and tensor current 



terms 5'!^'^\ Slj^^'^ and 5lj^^^ into a closed form expression we use the recursion relations 

i/j{n + l)=^{n) + -, V^(l) = -7^, ^'(n + l)=V;'(n)-^, ^'{l) = ^. (75) 



n 



n^ 



6 



Note that the three terms S^}'^', Slj;^' and 51}'^' stem from different diagrams. With the help of 
the results from Appendices A to C, resp., one obtains 



Si^ = A^ + 



9 






(76) 



The coefficients A^ and A^' are rational numbers. As an example we list the first values for A^ 
in the first parts of the three columns of Tabled] for rp = (scalar current), rp = 1 (vector 
current), and for the sake of completeness rp = 2 (tensor current, F = o^^\ We represent the 
moments in the form 






M 



N 



M] 



m 



N 



a 



1 + ^(61^-61 



IT 



Nj 



Mt 



< 



M, 



m 



M\ 



m 



N 



a. 



TC 



(K 






(77) 



As mentioned before, all moments defined in flTT]) are normalized to one in the sense that 
M^ = 1 for p*(z) = 1. Note that the difference 6n — Sn is scheme- independent. This feature 
was used in the high precision analysis of heavy quark properties [26] within NRQCD (see e.g. 
Ref. [2^). One can now easily find the actual magnitude of the correction. Indeed, for any 
desired precision and range of n, a set of perturbatively commensurate moments can be found. 
We therefore present differences of consecutive 5^ in the second part of the columns of Table [H 
Note that because of the s-integration in Eq. (!7T|) . moments represent massive vacuum 
bubbles, i.e. diagrams without external momenta with massive lines. These diagrams have been 
comprehensively analyzed in Refs. j2Hll^ . The analytical results for the first few moments at 
the three-loop level can be checked independently with existing computer programs (see e.g. 
Ref. BUI). 
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6 Applications to physics: Sum rules 

A phenomenological application is given by QCD sum rules (see e.g. Ref. [HI [32]) or finite 
energy sum rules (FESR, see e.g. Ref. [33]). In general, the main quantity of interest is the 
residue 

(|jB(0)|A(p,a)) = ABn(p,a) (78) 

where A{p,a)) is a baryon state and u{p,a) is the spinor that satisfies the free equation of 
motion (^ — mB)u{p, a) = 0. The contribution of the ground state baryon B to the spectral 
density of the correlator reads 

Pb(o)(s) = A^(5(s-m|) (79) 

while the excited states contribute to the spectral density Pb{s) that were calculated before. We 
define a threshold value sq where the excited states and the continuum start contributing. If we 
take moments, i.e. integrals over the phenomenological spectral density Pb{o){s) + 9{s ~ sq) pb{s) 
with different powers of s, we obtain the sum rule condition 

J^ PB{s)s''ds = J^ {pBm{s) + e{s-So)pB{s))s''ds = Xlm^j^ + J^ PB{s)s^ds (80) 

or, equivalently, the sum rule 

/■so 
>^Wb= PB{s)s^ds=:Mn{s^). (81) 

The two unknown parameters Sq and A^ {ttib is assumed to be known) can be determined in 
turn by the sum rule analysis. For this purpose we calculate the ratio of nearby moments, 

- - m^ (8^j 



7W„_i(so) /oVB(s)s"-irfs 
and adjust sq. Once sq is determined we can calculate A^ by using the zeroth moment, 

A| = Mo{so) = / pBis)ds. (83) 

Jo 

The second mass scale occuring in the problem is the mass m of the heavy constituent quark 
occuring in Pb{s)- It is fixed by assuming a value of 500 MeV for the difference A^ = tub — m 
between baryon and quark mass. Finally, we have to decide which of the sum rules we use. 
Here we have decided to use the n = sum rule. For p{s) we use Eq. ( ET\i in the scalar setting 
(A- type baryon: F = 1, F' = 75) or vector settings (E-type baryon: F = 7^, F' = 7^,^75, S*-type 
baryon: F = 7^^, F' = 1 with ^ — > '^^). ag is the running QCD coupling with a scale set by 
the heavy quark mass. 

6.1 A- type baryons 

The A-type baryon ground state masses are given by [1] 

m^+ = 2286.46 ± 0.14 MeV, m^o = 5624 ± 9 MeV. (84) 
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We first do the sum rule analysis for the LO contribution. Using the sum rule 



so 






/ ,P [s)(^s I Po{Z) — 



mi, 



(85) 



with zq = m^/so, we obtain the sum rule value for zq. Reinserting this value into J\4o{so), 
we obtain the value for A^+. The results are collected in Tab. O In going from LO to NLO, 



A^ 


zo (LO) 


Zo (NLO) 


A (LO) 


A (NLO) 


m-part 
g-part 


0.565550 
0.564191 


0.564193 
0.562752 


1.77978 X 10-^ GeV^ 
1.60523 X 10-=^ GeV^ 


3.01893 X 10-3 GeV^ 
2.72316 X 10-3 GeV3 




Ai; 


Zo (LO) 


Zo (NLO) 


A (LO) 


A (NLO) 


m-part 
g-part 


0.805403 
0.805088 


0.804903 
0.804575 


1.81465 X 10-3 GeV^ 
1.74339 X 10-3 GeV^ 


2.86564 X 10-3 GeV3 
2.75323 X 10-3 GeV3 



Table 2: Sum rule analysis for the A- type baryons with Ej\ 



rriA 



m 



500 MeV 



the value of zq does not change significantly. However, the values for A^+ and A^o change by 
a factor of 1.7 (A+) and 1.6 (A°) which disqualifies the sum rule method as an appropriate 
method to determine the baryon's parameters. Theoretically, we can improve the convergence 
behaviour by tuning the mass logarithm ln(?7i^//i^) to values of roughly 4.5 (m-part) or 6.0 
(g-part). However, the necessary values of /x (less than 200 MeV) are far from being realistic. 

6.2 E-type baryons 

For the Sg-type baryons, the three ground state masses read [1] 



mj.++ = 2454.02 ± 0.18 MeV, 



m^ 



2452.9 ± 0.4 MeV, m^o = 2453.76 ± 0.18 MeV. (86) 



We take the mass average m-£^ = 2.453.56 MeV as input parameter. For the recently found S^ 
state we take the average of the mass of the positive and negative charged state and obtain 
m-£^ = 5811.5 MeV [SB]. Tab. Ogives an overview over the results. The NLO/LO ratio for the 
residue is 1.6 for H^ and 1.46 for S^. 



6.3 E*-type baryons 

In order to extract the vector baryon components of highest spin 3/2 we use the Rarita- 
Schwinger formalism (see e.g. [37]). Instead of the spinor u{p,a) in Eq. (178!) we have to use 
the Rarita-Schwinger spinorM^(p, cr) (cr takes the four values ±3/2 and ±1/2). The equation 
of motion for the highest spin state leads to the two additional constraints 



7^M^(p,Cr) = 



and 



p^u^{p,a) 



0. 



^7) 
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Sc 


^0 (LO) 


^0 (NLO) 


A (LO) 


A (NLO) 


m-part 
g-part 


0.590168 
0.588942 


0.589363 
0.588086 


2.99916 X 10-^GeV^ 
2.72806 X 10-3 GeV^ 


4.73288 X 10-3 GeV3 
4.30472 X 10-3 GeV3 




s, 


20 (LO) 


zo (NLO) 


A (LO) 


A (NLO) 


m-part 
g-part 


0.811194 
0.810897 


0.810902 
0.810598 


3.10502 X 10-3 GeV^ 
2.98752 X 10-3 GeV^ 


4.53653 X 10-3 GeV3 
4.36492 X 10-3 GeV3 



Table 3: Sum rule analysis for the S-type baryons with E^, = mj] 



m 



500 MeV 



Using the ansatz 



n 



flU 



Aq^'q" + Bg""" + -a'^'^D^ 



al3 



where 

D^; = i?(<7^<7^ - gl^g';) + X,q^\glqp - g^q^) + X^g'^^'g^ - ^^g„), (89) 

we can contract from the left with 7^^ in order to satisfy the first constraint 7^11'^'^ = 0. We 
obtain Xi = -X2, A + AX2 = and 



R 



1 rA 



6 V2 



-Q 



B 



(90) 



Further imposing the condition g^II^'^ [38], we can determine all coefficients up to one for 
which we choose A. This gives rise to a basis of projectors (P^/^)^*^, ("P-j^^ )'^'', and (7^22 Y'^ ^^d 



nilpotent operators {V 



1/2- 
12 



liv 



flU 



(p3/2^) 

(^22 



l/2^/^^y 



^'P^'^^tMU 



12 



and (Vll^)''^ of the form ^ m\ 



fii/ 



9 

3 ' ' 



2 nf^n" 1 1 



3 g2 
IgV 



7^g'^) 



3 g2 



3?^^^^^ 



7V)i 



g- 
^ g"' 3g^ 



{V, 



l/2\/^i/ 



g2 3g' 



rg^(-g^ + 7' 



(91) 



In order to obtain the highest spin component we contract with [J^^/'^Y^ , If we use the above 
results for a sum rule analysis of the S*-type baryons with ground state masses P 



m^*++ 



2518.4 ±0.6 MeV, 



m, 



S•^ 



2517.5 ±2.3 MeV, m2*o = 2518.0 ± 0.5 MeV (92) 



with an average value of m^^ = 2518.0 MeV, we obtain the values in Table HI For the recently 
found S^-baryon we take the average of the mass of the positive and negative charged state 
and obtain m^- = 5832.7 MeV [36|). The ratio NLO/LO is 1.55 for S* and 1.45 for S*. 
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K 


zo (LO) 


Zo (NLO) 


A (LO) 


A (NLO) 


m-part 
g-part 


0.599077 
0.597895 


0.598274 
0.597045 


1.22636 X 10-^GeV^ 
1.11838 X 10-3 GeV^ 


1.90952 X 10-3 GeV3 
1.74130 X 10-3 GeV3 




n 


20 (LO) 


Zo (NLO) 


A (LO) 


A (NLO) 


m-part 
g-part 


0.811837 
0.811543 


0.811537 
0.811234 


1.26773 X 10-3 GeV^ 
1.21994 X 10-3 GeV3 


1.83958 X 10-3 GeV3 
1.77024 X 10-3 GeV3 



Table 4: Sum rule analysis for the S*-type baryons with E^,* 



mk 



m 



500 MeV 



7 Conclusions 



To conclude, we have computed the NLO perturbative corrections to the correlators of finite 
mass baryons containing one heavy quark and two massless quarks for a variety of quantum 
numbers of the baryonic currents. Technically this is a genuine three loop calculation with the 
two mass scales s and m?. We have considered both the massless limit and the threshold HQET 
limit of the correlator as special cases of the general finite mass formula. The two respective 
limiting expressions agree with previous massles and HQET results in the literature. From 
threshold to high energies the exact spectral density interpolates nicely between the leading 
order (leading in l/mq) HQET result close to threshold and the asymptotic mass zero result. 
This raises the hope that one can find a similar interpolation formula at the four loop NNLO 
level using the massless and the HQET four-loop results. These can be calculated using existing 
computational algorithms [lOl UT] . 
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A The light contributions 

The light contribution consists of contributions from the two diagrams (b21, light self energy) 
and (ell, light fish). These contributions are calculated in turn. 



A.l The light self energy (b21) 

The result has the same structure as the leading order contribution. 



Pb2l{S} 



16(47r)3^/' 



:^'^/'-'EpU(^V^)tr.(r,.) 



(Al) 



i=l 
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where trj(r, s) are given by Eqs. (!23l) and 



Pb2l[Z 



pLi(^ 



pLi(^ 



Pb2l{z 



PmA^ 



Ph2l{z 



e^ 



D-2 

2(D-1)2 

D-2 
2(D-1)2 

D-2 
2(D-1)2 

D-2 
2(1} -1)2 

L)-2 
2(D-1)2 

D-2 
2(1} -1)2 V£ 



1 



6£ 
1 

1 

1 

1 



17 
12 



Pal(^) 



3-7^+^)pri(^)-?f-+?^)prw 



p^i(^) 



12 
17 
12 



17\ .4*. ^ 



pri(^) 



17 
12 



|^)pS(.) 



2 /I 



slJ^I"-*^' 



^g* , 



(A2) 



A.2 The light fish (ell) 

For the hght fish the traces are more comphcated. However, using the parameter cr introduced 
in Eg. (1361) . the result again has again the same structure as the leading order contribution, 



where 

Pcll(^ 



Pell (2; 



^cll^ 



Pcll(^ 



p: 



cll^ 



Pell 12: 



PclllS) 



9lG^ ,.3^/2-^EPcii(^V^)tr.(r,.) 



16(47r)3^/2 



i=l 



^^ + 19-24C(3)4(A + ^-l 



^19-24C(3)4(A + ^_1 



- + 19-24C(3) + ^f-| + --l 



12 \e 



+ 19 - 24C(3) + 



12 \e 



%?-i 



j + 19-24C(3)4(| 



6 



- 1 



^ + 19-24C(3)4(A + ^-l 



Pal(^), 



Pal(^), 



P^l(^), 



P^l(^) 



2 + ci 



3£ 



2 + cM^ + l 



pr(^), 



(A3) 



pri(^). (A4) 



Note that the parameter cr only appears in squared form, so that the signature of the current 
(i.e. whether or not there is a 75) is irrelevant for the result. 

A. 3 Merger with the leading order diagram 

When one combines the results for the two light diagrams and adds the leading order contri- 
bution, the divergences can be factored out. When one combines the contributions one has to 
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take into account colour and combinatorial factors. The relevant factors are A''^! for the leading 
order diagram (al), 2Nc\Cf for the light self energy diagram (b21), and —NJCb for the light 
fish (ell). The results read 

/ \ TIT 1 / \ TIT 1 [^ — 2)LrS 1 ^ , 2 / \ 

Pleading (S) = Nc\pal{s) = ^'^- IQU^\D fJJ _ 1)2 "P^l ("^ /*)' 






PHght(^) = -^ir^\(-j + - + B2]plr{s)+[^ + B',]p':,{s) 



where the coefficients Bi, B'- are given in Eq. (1351) . and 

6 

pai{m^/s) := ^pli(mVs)tri(r,s), 
1=1 

Pliim'/s) := 5:pri(mVs)tr,(r,5), 

p':,{mys) := p:{%mys)tT,{T,s)+plUmys)MT,s). (A6) 

The functions Pai{z), p*i(z), and Pai{z) that occur in the spectral densities for the leading and 
light contributions vanish for e —>■ 0. This can be seen by retracing the construction down to the 
elements Pai{z) and p"i{z) {i = 1, 2, . . . , 6). These elements are given by a linear combination 
of spectral functions pv{l,ne — p; z) as defined in Eq. ( !29l) . Using the abbreviations 

^ ^ ^ 1 ^p ^ T{ne-l)T{3-{n + l)e) 

" r(n£-l)r(3-(n + l)£)' ^' r(n£-p)r(p + 2-(n + l)£) ^ ^ 

for the overall and relative factor, the elements can now be written as 

pIi(^) = Cisgiiz), pi\{z) = C2egi{z), 

p'Uz) = C,eg'l{z), p':iiz) = C^eg'^iz), ^ = 1, 2, . . . , 6 (A8) 

where (note that C^ = 1) 



9n{^) = DLlgn2{z) + C',{l-z)'gno{z)-2(^^^~zyg^,iz)V 
9ni^) = -|c^^n2(^)+C°(l-^f^„o(;^) + 2(^^ + z)^„i(;.)|, 
g'niz) = ^{clgns{z) + C'^[^ + zyi-zrg^o{z) + 



^^ /6-D \ , , [2- D 2Dz ,, , 
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g'niz) = -Uc^gnsiz) - C°(l - zfgnoiz) + Cl{l + ?,z)gn2{z) - (1 - z){l + 3z)^„i(2;)|, 
gl{z) = -\\clgn^{z) - C°(l - zfgr,^{z) + Cl{l + ^z)gr,2{z) - (1 - 2;)(1 + 32;)^„i(^)|, 

9l{^) = -ll^Clg^3{z) + C'^{l-zni + z)g^o{z) + 

^2/7D-10 ,\. ,, /7D-10 ,3D-4 .\ ,,1 

+ " [ D-2 + j ^"'^''^ "^ ( ^_2 + ^^D^ry^ + 3^ j ^ni(^)|. 

^;'(^) = 9n{^) = 9n{^) = g'niz) = 0, 
9n{z) = -9nl{z), 






9n{Z 



U{l + z)gn,{z)+Clgn2{z)}. (A9) 



We can cast this into a more compact form by writing 

pS {z) = CieM'^ (z) , pit (z) = C^ecf^ (z) ( AlO) 

where 

g^Krn'/s) = ^^f (mVs)tr,(r,s). (All) 

These expressions can then be inserted into the sum of leading and light contributions. One 
obtains 

GN,\ 2-2. (^ - 2) ^ 

Pl{s) := Plcading(s)+ Plight (S) = J^T^^^ (0-1^ ^ ^ 



' 4. ^'^ij(Z^3%^(^^ + ^^^)^^(-V.)| - (A12) 



16(47r)^ (D-iy 'Y'' ' ' Ane ^'^ ' ' ' 47r 

The coefficients C" and C"^' can be expanded in e, viz. 

C"- := s--^(5o + 5i6 + 52£2) = C"-(°)+C""(i)£ + 0(£2), 

C- = ^-^^^^^^(^1 + ^2^) = C7-(°) + 0(.). (A13) 

The main singularity of the NLO contribution is proportional to ^2 (-2)- But ^2 (-2) is similar to 
^i(z). Indeed, one obtains 

hiiz) := Uz) - h{z) = gi\\z)e + 0{e'). (A14) 
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Therefore, the leading singular NLO part is added and subtracted where §2(2) is replaced by 
gi{z). The leading singular part of this new contribution is then combined with the LO part. 
Expanding 

g^iz) = gf\z) + 0(5), g',{z) = gT\z) + 0{e) (A15) 

and finally using D = 4, one obtains 



GN,\ 2-2e (£-2) r / asCi. (o)\ . , 2/ ^ ^ 






16(47r)^ {D - 1) 

9(47r)4 (V 47r£ / ^^ ' ' 

+ ^ (C^^^^)#K/^) + C^^'^g^hm'/s) + C-(°)^r(mV.)) }. (A16) 

The LO contributions p\i{z) are given by Eqs. fIST]) . For the remaining light contributions one 
has 

pU(^) = ^(2 - 4rr + r^)/3U^) + ^^rPll(^) , (A17) 

where 

^w X 19 41 29 2 13 . 29 4 / 10 . 2 o s 1 4A , /. 

p;^(2) = — + —2 z^ + —z^ z^ - [l + —z-Qz^ + 2z^ - -z^]\Yi{l- z) + 

'^''^ ^ 24 18 6 6 72 V 3 3 / ^ ^ 

+ Q + ^2 - 62' + 2^3 - }-z'\ hiz + Az {hi2{z) - Li2(l) + ^ In^ z") , 

.9/ X 25 41 47 2 ^ 29 4 /I 11 2 s 1 4\ , / 

pI(z) = — + —Z z^ - z^ + z^ -[- + —z-^z^ - z^ + —z*]\n{l- z) + 

'^''^ ' 96 9 12 288 V4 3 12 / ^ ' 

+ Q + y^-^^'-^' + ^^')ln^ + (2^ + 3^2)(Li2(^)-Li2(l) + iln2;.), 

.0, , 107 59 47 2 67 o 149 . 157 , 

piiz) = z-\ z^ z-^ -\ z^ z^ + 

^"^ ' 600 120 30 30 120 600 



^ ~ ,2 o~3 



-[--z + 2z' - 2z^ + z^- \^z^\ ln(l - z) + 

+ (— - -^ + 2^2 - 2z^ + z^- -z^") Inz, 

^4/ N ^5/ N ^6/ ^ 137 107 47 2 9 3 47 4 157 , 

p'^iz) = pI{z) = p^(z) = z ^2 H z^ 2^H z^ + 

^"^ ^ ^"^ ' ^"^ ' 2400 240 180 10 160 3600 

z z^ + z^ z^ H z^] ln(l - z) + 

20 2 3 4 30 / ^ ^ 

+ ( z z^ + z^ z^ H z^] hiz- z^ i^hiz) - Lisfl) + - In^ z 

\40 4 3 430/ 1^2v; 2v;^2 



(A18) 



29 



B The self energy correction of the massive line (bll) 

The first order self energy correction of tlie massive line is given by the inclusion of the master 
bubble 

^B{k') = I fJnn. L n. = ^?^I^V^(1, 1^ ^V-^^) • (Bl) 



. (27r)^(/2-m2)(A;-/)2 {Att)^/^ 

The corresponding spectral density reads 

/^2\D/2~2 5-^^/2-2 

Pnis) = ——^j^pvil, 1; s/m^) = —--^pv{l, 1; m^/s). (B2) 

With the help of the dispersive representation it is easy to see that the mass and momentum 
part of the self energy correction can be written as 



S„(A;2) = Dg^.J 



PB{s)ds 



2\ 2 — -D 2 /" ( 1 , m-] PBis)ds 



In order to enact the renormalization and to absorb the singular parts of these NLO contribu- 
tions in the renormalization factors for the mass and the wave function, we use an expansion 
for the propagator-type factor, 

- — ^ — - — —r-^-r, ^ T, [ -ika-imih -a)]- . (B4) 

One obtains 

2\ ^ n2\ 2-/} 2 t PB{s)ds ( m2\ r pa{s)ds 






s-A;2\^2 2 s J ' s-k^ 

Using momentum subtraction at k'^ = m'^, the singular parts can be split off, 

.j^2. ^ f Pa{s)ds ^ f Pa{s)ds M 1 ^_\ 

J s — k"^ J s — m? J \s ~ k"^ s^ — m? J 

Pa{s)ds 2 2. [ Pa{s)ds . / 2n , /,2 



s — rn? J [s — w?){s — k"^) 



where 



2a / Pa{s)ds -gl f fi'^Y G 
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Having absorbed the divergent parts into the renormahzation of mass and wave function, the 
finite parts can be expanded again, 

i(l + af) iJjt ( 2m%f \ im ( 2m%f \ 



^ — m{l + hf) k"^ — rn? \ k"^ — m? ) k"^ — m? \ / / )^2 _ ^ j 

'^ ^1 + P(A:^)) + T,^ (l + M(A;^)) . (B8) 



For the leading order diagram one obtains 

^, , 2n 2^(1,1) /• d'^k i ( , ,,2nD/2-i ^^ ^ 

where (7(1, 1) = Gje is the massless master bubble. The corrections to the propagator results 
in the use of an effective propagator 

«-(*') = f3^««(*^) + p^«o"^(*') (BIO) 



with (pa+fe(s) = Pa{s) + Pfe(s)) 



ifc (u2^ _ ' • / Pa{s)ds 2 /■ Ph{s)ds ( 1 1 



® k'^ — vn?' J (s — m?){k'^ — s) J {s — m?Y \^^ ^ ""^^ k^ — sj ^ 

,^,,2x ^ , /■ Pa+bis)ds 2 /■ Pfe(s)c?S / 1 1 



k'^ — m^ J [s — m'^){k'^ — s) J (s — -m^)"^ V/e^ — m^ k'^ — sj 

One finally obtains 

Ki(g ) + Hn(g ) - -^^^^pT^ j (^2n)^ ^-{g - kyy-^/^' ^^^^^ 

It is therefore obvious how to calculate the contribution Vbii{q'^) from the self energy correction 
of the massive line and its spectral density pbn{s) = ^Pliii^) + "^Pwil-^)- After integrating by 
parts the final result can be seen to be a convolution of the leading order contribution with 
specified weight functions 



where 

^"*'^'* =• /4^N'z)/2 P"("^Vg). Pa+6(s) =: 7^^^f^Pa+6(mVs) (B14) 

and 

^'(^) = f (?^''^' =^ (i|vji'('"V.). (B15) 
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The corresponding spectral functions are 

„i , , 49 116 74 2 28 3 59 4 /I 44 ,24324. 

nl-.Jz) = 1 z z^ H z'^ z^ + - H z - Az^ + -z'^ z^ In z + 

mn ; 36 27 9 9 108 V3 9 3 9 ' 



2 20 
3 



+ ^z - Az' + ^z' - ^z'^ ln(l -z) + 2z ^U,{z) - U,{1) + ^ In^ z^ 



.0 , , 55 329 46 2 3 3 59 . / 1 34 3 2 2, 1 , , , 

HbuK J 144^54 9 2 432 V12 9 2 3 IP ' 

1 22 2 2 3 1 A , , 

- + —z~2z^- -z^ + —z^ ln(l -z) + 
6 9 3 18 / ^ ^ 

+ ^z(2 + 3^) (^Li2(z) - Li2(l) + ^ In' ^ 

.3 , , 77 433 289 o 299 3 613 . 623 . 

pf (2) = z-\ z^ z-^ -\ z* z^ + 

^^^^^ ' 1800 360 90 90 360 1800 

f ^ 5 42 43 24 2 5\ 

V -z z^ + -z-^ z^ H z^ In ^ + 

V30 6 3 3 3 15 / 

/2 2 42 43 24 2 5\,^ , 

z + -z^ z^ + -z^ z^ ln(l - z), 

V15 3 3 3 3 25 / ^ '' 

.4 , , .5 , , .6 , , 107 709 173 2 91 3 199 4 623 5 

ptAz) = PmAz) = PMTiz) = z z^ z'^ z^ ^ z^ + 

miv ; mu ; miK > -^200 720 1080 60 480 10800 

(I 5 47 2 23 z^ z^\^ 

1 z H z^ z'^ ^ In z + 

V120 12 36 3 6 45/ 

( ^ ^ 2 2 2 3 z^ z5\ , 

z^ + -z-^ \ n 1 - z) + 

V30 3 9 3 6 45/ ^ ' 

-^^2(Li2(^)-Li2(l) + iln2;.). (B16) 

These results have already been renormalized. In addition to Eq. flB13P we have to take into 
account a further finite contribution coming from the singular parts of a and h. Since 

_ 2 / ,,2\^ n . „2 



^)'fO-^0(^^')' - + ''-(4^(^)'f(2-2^ + 0(^^)) («1^) 



' (47r)^/2 \m 

we can absorb the singularity into the renormalization factor. This is the case for the MS-mass 
where the finite constants have to be added to the result. If we absorb the finite constants as 
well we end up with the pole mass. This is preferable because then we do not have to take care 
of the numerator singularity containing h. In any case, we obtain expressions for ln(yU^/m^) 
which have the same coefficients as the poles. 

C The semi-massive fish contribution (c21) 

In order to determine the semi-massive fish contribution one has to calculate a number of scalar 
two-loop integrals. These so-called prototypes are spectral functions pv{ni, 77-2, ^3, 77,4, ns; s/rn?). 
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Figure 3: Different representations of the tetrahedron or proper fish topology (a-b) and the 
spectacle topology (c-d) 



These spectral functions are determined by the discontinuities of the correlator functions 

V{ni, n2, ns, Ua, n^, q^/rn^) given by 



1 



(4vr) 



D 



,ra 



2\D— ni— n2— ng— ?i4— ns 



V"(ni, 722, ng, 724, Tig; q^ / m^) 



d^'k d^^l 



Dl 



(27r)^ (27r)^ (P + m'^Y^{P + m^Y^{{q - kfY-^{{q - lYY^{{k - IfY^ 



(CI) 



(for convenience written in the Euclidean domain). For later use it is convenient to use the 
representation 



Py(ni,n2,n3,n4,n5;zj 



z 



D—ni—n2—n3—n4—n^ 



Py(ni,n2,n3,n4,n5;l/z). 



(C2) 



A subset of the prototypes turn out to be reducible to scalar one-loop integrals using the 
spectral representation pv{ni,n2; z) := z^/'^~'^'^~"'^pv{ni,n2;l/z), where the spectral function 
Pv{ni,n2', s/m?) is given by the discontinuity of the correlator function V{ni,n2'i q^ /m^). One 
obtains 



[m 



2\D/2-ni-n2 



V{ni,n2]q^/inn? 



d'^k 



1 



(47r)^/2' ' V 1, .,^ / / • y (27r)^(P + m2)"i((g-A;)2)"2 • 

An alternative approach is to relate them to the massless one-loop integrals G{ni,n2) with 



(C3) 



(47r)^/- 



■(? 



2\D/2-ni-n2 



G{ni,n2) 



d^k 



(27r)^(F)"i((g-A;)2)"2 ' 



(C4) 



C.l The proper fish and the spectacle prototype 



We start with the most difficult prototype py(l, 1, 1, 1, 1; z), the proper fish or tetrahedron (see 
Fig. [3](a) and (b) for two different representations of this topology). The corresponding correla- 
tor function can be obtained as a limiting case of an expression taken from the literature [lllll2]. 
The discontinuity of this expression turns out to be finite. One obtains 



(C5) 



pvil, 1, 1, 1, 1; z) = 4 ( Li2(z) + - ln(l - z)\nz] + 0{e). 
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If the last entry "1" is replaced by "0", the correlator consists of two master bubbles. This 
diagram is termed the spectacle diagram (cif. Figl3](c) and (d)) where the name derives from 
the pictorial representation Fig 131(d). We calculate the spectacle diagram by adding a further 
scalar line, 

J^ni, 1,1,1, 0, 1; qVm^) = J (0J (^ / ^(1' 1; pV-^)^(1, 1; pVm^ (C6) 

Note, however, that we have to subtract the infrared divergence before we can make use of the 
dispersive representation. We insert 

V{l,l;p m)= ds = V[l,l]-l) - [p + m ) C7 

J s + p'^ J [s — m^)[s + p'^) 

where we have chosen the momentum subtraction at the point p^ = —m? in the Minkowskian 
domain (cf. Appendix B). One obtains 



\ni 



2U-e 



^a(g^ si)pv{l, 1; si/m^)dsi + / Afe(g^ si, S2)py(l, 1; si/m^)pv{l, 1; S2/m^)dsids2. 

The convolution functions are given by 

n f d^P V^(l,l;-l)(p2 + m2) ^^^^^ ^^ ^ ~ ^ 2 



Aa(g,si) = -2/——^- 4^- -^— = 2^(1, l;-l)Ao(g ,si). 



A6(g^Sl,S2) = j 



d^p {p^ + m^) 



2\2 



(27r)-^ (si — iv?){s2 — ni?){si + p^){s2 + p'^){q — pY 

1 r d^p ( ijP' + m?Y (jp' + m^Y 



(C9) 



{si - iv?){s2 - ■m'^){s2 - Si) J [271)^ \{si + p^){q - pY {s2 + P^){q - pY , 
{si - m^)^Ao(g^ si) - {.$2 - m^YMQ^, ^2) 

(Si - W?){S2 - 'W?){S2 - Si) 

where 

Note that for the reduction to the fundamental convolution function Ao((3'^, s) one can use the 
fact that each non-negative integer power of p^ occuring in the integrand of this function can 
be effectively replaced by —s. Calculating the spectral function we obtain 

|^^pv(l,l,l,l,0,l;s/m2) = (Cll) 

= / Aa(s, 51)^^(1, 1; si/m^)dsi + / Xb{s, si, S2)py(l, 1; si/m^)py(l, 1; S2/rrY)dsids2 
where 

Aa(s, si) = 2\/(l,l;-l)Ao(s, Si), Ao(s, si) = — — ^py(l, 1; s/si), 

Afe(s,si,S2) = ^ j-7 j-7 7 . (C12) 
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In order to combine both parts into an integral including a unique convolution function and a 
single integrand (which later on will be identified with the prototype pv{^, 1, 1, 1, 0; s/w?)) we 
use the fact that the second integral is symmetric in si and S2- We further make use of the 
expression 

at D = 4 (i.e. e = 0) for the one-loop spectral function to perform one of the integrations as a 
principal value integral. One obtains 



Afe(s,Si) := I \b{s,Si,S2)pv{lA]S2/m^)ds2 = 

= -2A„(s..0^i^^ln(ii^). (C14) 

such that 

^^-^pv(l,l,l,l,0,l;s/m2) = 

2Ao(s, si) (v{l, 1; -1) -(l-^] In (^ " l)] Pvih 1; s,/m')ds,. (C15) 



We identify 



pv{l, 1, 1, 1, 0; l/z) = 2 (V{1, 1; -1) - (1 - z) In Q " l)) Pv{l, 1; 1 A) (C16) 

and finally obtain 

pv{l, 1, 1, 1, 0; z) = 2\/(l, 1; -l)z-'pv{l, 1; z) + 2(1 - zf (ln(l - ^) - In^) . (C17) 

Next we calculate the prototype py(l, 1, 1, 1, —1; z). For such prototypes with negative entries 
we need the vector integral V'{1, l;p'^/m^) defined by 



(47r)^/2' ' /^ V > ./^ / / y (27r)«(P + m2)(p-A;) 

One obtains 

V'HA-v'/m') = - ( 1 - — j \/(l,l;pVm2) - _K(1,0;-1). (C19) 

Again, we use momentum subtraction at the point p^ = —m'^, 



V'{1, l-p'/m^) = \/'(l, 1; -1) - (p2 + m^) / ^^'.'"^ %: (C20) 

J (s+p^j(s — m^j 

where the spectral function is given by 

1 / m'^\ 1 

p'yil, 1; s/m2) = - h + _ j p^(l, 1; s/m2), p'^(l, 1; ^) = -(1 + ^)p^(l, 1; ^). (C21) 
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In terms of the above vector integral we calculate 



V(l, 1, 1, 1, -l-y/m') = -2\/(l, 1, 1, 1, 0;pVm2) - 24v^'(l, 1;pV^')', (C22) 



and obtain 

py(l, 1,1,1, -1; z) = -2zpv{l, 1, 1, 1, 0; z) + 4V'{1, 1; -1)^-%(1, 1; z) + 

-(l + z)(l-^)2(^l + (l + z)ln(^i-l)) = 
= -W{1, 1; -l)^i-^/)y(l, 1; ^) + 4\/'(l, 1; -1)^-%(1, 1; z) + 

-(1 +z){l- zf - (1 - zf (ln(l - z) - Inz) . (C23) 

C.2 Prototypes of the class /3y(l, 1, 0, 1, 1; z) 

Prototypes with one vanishing massless propagator reduce to a nested integral. For the general 
case that we need to consider here we obtain 

1 [^^^)^-^-n.-~n,y{l^l^^^n,,n^) = 



(47r)^ 

d^k 1 r d^l 



V{\, n^; k /m ) = V{1, n^;-!) - [k +m ) — — -- (C25) 



(27r)^ (fc2 + m^){{p - A;)2)"4 J (27r)« {P + m^){{k - /)2)"5 

= 1 /^2^D/2-l^n. f ^^A; V (1, 725; ^Vm^) 

Again making use of momentum subtraction 

(s — m?){s + A;2) 
one ends up with 

V^(l,l,0,n4,n5;pV^2) = V{l,n,-,-l)V{l,n,-y /m') - j PX^}^J^^hlhlflv{l^n,-y / s)ds, 

J s — m'^ 

Pv'(l,l,0,n4,n5;s/m^) = V{l,n^]-l)pv{l,ni]s/m^) - — — \ V{l,n4; s/si)dsi 

J s — m^ 

and (C26) 

pvil,l, 0, 714,7x5; z) = 

Because of py(l, 0; z) = we obtain 

py(l,l,0,l,l;z) = \/(l,l;-l)^-^py(l,l;^) + l-z + lnz, 

pv{l,l, 0,1,0; z) = V{l,0;-l)z'''pv{l,l;z) + {l-z)z, 

pv{hl,0,0,l;z) = (C28) 
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(furthermore, pi/(l, 1,1,0, ns; z) = pvi^, 1, 0, 1, n^; z) because of the symmetry of the problem). 
For the last prototype of this class, /^^(l, 1, — 1, 1, 1; z), the vector integral V will appear again. 
Using a dispersion relation we obtain 

py{l, 1, -1, 1, 1; 2;) + (1 + z)py{l, 1, 0, 1, 1; z) = 

= 2z'^V\l, 1; -1)pU1, l;z)-2 [' ^Hilli^p'^(l, 1; z/z,)dz,. (C29) 

J z Zi[i — Zi) 

The final result reads 

pv{l, 1, -1, 1, 1; ^) = (1 + z)V{l, 1; -l)z-'pv{l, 1; z) + 2ni, 1; -lK^Pv(l, 1; ^) + 

-l + z+\z'-(^l + zynz. (C30) 

C.3 Prototypes of the class /3y(0, 1, 1, 1, 1; z) 

If one of the massive propagators vanishes, the result is given by the product of a massive and 
a massless one-loop correlator. The spectral function reads 

py(0, n2, ng, 724, u^, z) = G{n3, n5)py(n2, ng + n4 + n5 - L'/2; 2;). (C31) 

For the special cases that occur in our calculations we obtain 

pv(0,l,l,l,l;z) = G'(l,l)pv(l,£ + l;z), 

pv(0,l,l,0,l;z) = G'(l,l)pv(l,5;z), 

pv(0,l,l,-l,l;z) = ^(1,1^(1,6-1;^), (C32) 

as well as pv{0,0,n3,n4,n5; z) = 0. The last prototype pv'(— 1, 1, 1, 0, 1; z) is more difficult. 
With a little bit of work one finds 

py(-l, 1, 1, 0, 1; z) = -^G{1, 1) ((1 - z)pv{l, e; z) + pv.(l, e-l;z)). (C33) 

In order to calculate the final result, we consider the remaining spectral functions 

P.(M-M) - r(.-i)r(3-2.) r(^-^""-^"(--)""-^-' 

Pv{l.e-z) = r(e)r(2 - 2e) l! ^^ ~ ^)'''^^^''(^ " ^Y^'^d^^ 
py{he+l;z) = ^^^ ^ ^^^^^ _ ^^^ J\l - xr'^x%x - z)-'^dx. (C34) 

The first two integrals can be evaluated for e = 0, while for the last member of this family the 
singularity in G(l, 1) is not cancelled. However, we can subtract and add pv{^, 1; z) to separate 
the singular and finite parts. Using 

py(l, 1; z) = YrTTe) j^^ ~ "^^''^'^ " ""^"''^'^ ^^^^^ 
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and r(l + e)T{l - 2e) = T{1 - e) + 0{e^), we obtain 



-[^pv{l,e + l;z)- z pv{l,l\z, 
e Jz 



(1 - xY'{x - z)-' ((1 - x)-^x^{x - zy - z"^) dx + 0{e) 
/ (In z — ln(l — x) + In X — ln(a; — z)) dx + 0{e) = 
l-z + {l-2z)lnz-2{l- z) ln(l - z) + 0(e). 



(C36) 



C.4 Table containing all needed prototypes 

All necessary prototypes are listed in this subsection, starting from the most complicated one, 
the proper fish prototype, to those that are zero. Using 0(1, 1) = Gje and 



1^(1,1 
1^(1,0 



-1) 
-1) 

-1) 



Vie) 



G 



\-2e 
V(e-\ 



+ 0(e) 






(C37) 



one has (in addition to inherent symmetries) 



py(l,l,l,l,l;z) = 4(^Li2(^) + -ln(l-^)lnzj, 
py(l,l,l,l,0;z) = 2^z-^py(l,l;^)-2(l-^)2(ln(l-^)-lnz), 

py(l, 1, 1, 1, -1; z) = (1 - 3z)-z-'pv{l, 1; z) + {l- z^)z - (1 - zf (ln(l - z) - \nz) 

G 

-- -jz^-'pv{l,l;z) + {l-z)z, 



py(l,l,l,0,0;z) 
Py(l, 1,-1,1,1;^) 



--(1 + 2;)-2; ^pv(l, l;z)-- + z--z^-(^- + z]lnz, 



G 



py(l,l,0,l,l;z) = —z''pv{l,l',z) + l-z + \nz, 

py(0,l,l, l,l;z) = -z-'pvil,l;z) + l-z + {l-2z)lnz-2{l-z)ln{l-z), 

I- z^ 
py(0, 1, 1,0, l;z) = — \-z\nz, 



2 

-(- 

2 V3 

1 1 2 

\- z r 

6 2 



3^ - 3z' 



-z^ + 2z{l + z)\nzj , 



1 



2;^ + z^\nz 



py(0,l, 1,-1,1; 2;) = 

py(-l,l,l,0,l;z) = 

while 

py(l,l,0,0,l;z) = py(0,l,l,l,0;z) = py(0, 1, 0, 1, 1; z) = py(0, 0, 1, 1, 1; z) = 



(C38) 
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pv(l,l,0,0,0;2;) = /)y(0, 1, 1, 0, 0; z) = py(0, 0, 1, 1, 0; 2;) = pv(0, 0, 0, 1, 1; 2;) = 0, 

pv(0,l,0,0,l;^) = py(0,l,0,l,0;z) = 0. (C39) 
The symmetries are given by 

py(n2, ni, n4, ng, ng; z) = pv(«i, ^^2, %, ^4, n^, z). (C40) 

C.5 Spectral functions for the semi- massive fish 

The spectral functions that we have obtained in the course of our calculation are 

.1 , , 49 67 569 2 169 o 25 . /I 14 43 2 32 . 25 A, 

pi (2) = \ z z^ H z^ z^ + - H z z^ H z-^ z^ ]\nz + 

'^"^^^ ' 36 54 108 54 54 V3 9 9 9 54 7 

/31 26 . 2 34 3 25 A, ,^ 

\ z-^z^ A / / ln(l - 2) + 

V18 27 9 54 y ^ ^ 

+ g + ^.-8.^ + ^.^-^.^)(Li,(.) + iln.ln(l-.)) + 

^iUz-z^^ ^) (Li2(^) - Li2(l) + \ In^ ^) + 8z (^Li3(z) - 113(1) - 1 lnzLi2(z)) , 

.1/ . X 25 13 49 2 17 3 -2^ /I 5 2 8 3 z^\ , 

^-(^) = -144-54^+36^ -18^ -432 -(12 + 9^-2^ +r+¥J^^^ + 

+ (^ + ^^ - 2.^ + ^.^ + ^) ln(l - .) - ^.(1 - Z-) (Li2(.) - Li2(l) + \ In^ .) , 

. , , 55 665 895 2 325 3 25 . 

p^ (2) = \ z z^ r H z^ + 

Pc2iv ; ^44 216 432 216 216 

( \ 35 77 2 11 s 25 A , 

+ \ z z^ z-^ H z^ In 2 + 

V12 18 36 6 216 / 

37 115 ^2 17 3 25 A, ,^ 

\ z - ?,z^ z^ H z^ n(l - z) + 

72 27 9 216 / ^ ^ 

+ (^^ + ^z^ + ^^' - ^) (Li2(^) - Li2(-1) + ^ In' ^ 
+ {Az + 6^2) (^Li3(^) - Li3(l) - i ln^Li2(^; 

.2/ / X 31 41 5 2-2^ -2^ /I 11 -2^ 7 3 ^^X , 

'^"'^^ ' 576 54 16 2 1728 \AS 18 8 9 36/ 

+ [^Q + l'- l^' - 1^) 1^(1 - ^) - I (1 + 3^ + ^') (Li2(^) - Li2(l) + i In' ^ ) , 

n , , 257 349 1273 1321 3 1507 . 58 , 

p;^ [z] = z H z^ z"^ H z^ z^ + 

"^"^^^ ' 900 900 900 300 450 225 
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/I 2 41 2 27 3 101 4 17 5\, 

V15 5 15 5 30 450 / 

/13 3 34 2 ^ 3 7 4 17 5\, ,, 

2 H z^ - 6/ + -/ z^ n(l - z) + 

V50 2 9 2 450 / ^ ' 

+ (-- -z+-z' - -z' + -z* - -z'>) (Lhiz) + -lnzln(l-z)] + 
Vl5 3 3 3 3 15 / V ^ ^ 2 ^ V 

+ (-2^^ + A^^j (^u,{z) - Li2(l) + i In^ z) , 

.,, , , 107 77 107 2 7 o 347 4 7 . 

p;^' (^) = \ z z^ H z-^ z* H z^ + 

' ''^^^ ' 3600 720 360 10 720 3600 

I \ z z^ 23 13 4 ^n , 

\ z-^ A z* H In 2 + 

V60 12 3 3 36 15; 

. , , . , , 287 124 16871 2 3097 , 1639 4 29 , 

p^ (z) = p% iz) = z z^ H z'^ z^ H z^ + 

Hc2i\ I Hc2i\ ) 3gQQ 225 10800 1350 5400 675 

(I 2 217 2 88 z^ 17 A , 

+ ieO - I5= - I80^ +i5^ +120 + 2700= j'"^ + 

/ 97 31 35 2 19 , z^ 17 A , ,^ 

— z z^ H z'^ \ z^ ln(l - z) + 



V1800 36 27 9 72 2700 

( ^ 2 4 2 43 z^ 2 5\ / . , , 1 , ,^ 

+ 2 z^ + -z-^ \ z^ Lisf^) + -\nz 11(1 - z] 

V15 3 9 3 3 45 ; V ^ ^ 2 ^ 

+ (4--l-4-g(L.W-L.,l)4ln^. 

-2z^ {hh{z) -hh{l) -UnzU^^z] 

-ii ( \ - -5' f \ - -^6/ / N _ 137 131 _ ^ _ 11^ 3 1^ 4 ^ 5 

Pc2i(^j - Pc2i{^) - Pc2i{^) - ^44QQ + ^44Q^ 4320 730"^ + 960^^ 21600^^ + 

/ 1 z 7 ^ z^ 19 , z^\, 

-z + — - — -7-2 - 777: In -2 + 



V240 24 72 3 144 90^ 

/I 2 5 2-2^ 19 4 -2^\w. 
+ 2^ H z^ ln(l - z) 

V240 18 36 3 144 90/ ^ ^ 

+ ^{2-z')(Li,{z)-U,{l) + hn'z 

. , , 287 221 3821 . 3569 , 1141 4 29 , 

^0 (^r) = z z H z z H z + 

'^'^^^^ ' 3600 900 10800 2700 1350 675 
/I 2 53 2 98 3 367 4 17 A , 

V60 15 180 45 360 2700 J 
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/49 13 8 2 7 3 25 4 17 .X , ,, 

z z^ + -z^ z^ H z^ ln(l - ^) + 

V600 12 27 3 24 2700 / ^ ' 

/I 2 4 2 43 z^ 2 gX / . , , 1 , ,^ A 

+ (4^4-g)(L.2(.)-L.(l) + iln^.) + 

- 2z=^ (^x^{z) - Li3(l) - i ln^Li2(^)) . (C41) 
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